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Abstract 

Because of flux conservation, flux- averaging justifies the use of the distance- 
redshift relation for a smooth universe in the analysis of type la supernova (SN 
la) data. We have combined the SN la data from the High- 2; SN Search and the 
Supernova Cosmology Project, and binned the combined data by flux-averaging 
in redshift intervals of Az = 0.05 and Az = 0.1. We find that the unbinned data 
yield a Hubble constant of ifo = 65 ± 1 kms^^ Mpc~^ (statistical error only), a 
matter density fraction of Qm = 0.7 ± 0.4, and a vacuum energy density fraction 
of Oa = 1.2 ±0.5. The binned data for Az = 0.1 yield Ho = 65 ± 1 kms"^ Mpc"^ 
(statistical error only), flm = 0.3 ± 0.6, and = 0.7 ± 0.7. Our results are not 
sensitive to the redshift bin size. Flux-averaging leads to less biased estimates of 
the cosmological parameters by reducing the bias due to systematic effects such 
as weak lensing. 

Comparison of the data of 18 SNe la published by both groups yields a mean 
SN la peak absolute magnitude of Mb = —19.33 ± 0.25. The internal dispersion 
of each data set is about 0.20 magnitude in the calibrated SN la peak absolute 
magnitudes. The difference in analysis techniques introduces an additional 
uncertainty of about 0.15 magnitude. 

If the SNe la peak absolute luminosity changes with redshift due to 
evolution, our ability to measure the cosmological parameters from SN la data 
will be significantly diminished. Assuming power-law evolution in the peak 
absolute luminosity, (1 + z)^, we find a strong degeneracy between the evolution 
power-law index j3 and the matter density fraction Q^n- For = 0.3, we 
find that the unbinned data yields i^o = 65 ± lkms~^ Mpc~^ (statistical error 
only), ^A = 1-4 ± 1.1, and f3 = 0.5 ± 1.6, and the binned data (with Az = 0.1) 
yields Hq = 65 ± lkms~^Mpc^^ (statistical error only), = 0.6 ± 1.4, and 
/3 = 0.0 ± 1.0. 



^Present address: Dept. of Physics, 225 Nieuwland Science Hall, University of Notre 
Dame, Notre Dame, IN 46556-5670. email: Yun.Wang.92@nd.edu 
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1. Introduction 

Near the end of the millennium, cosmology has become a data-driven science. We are 
closer than ever to determining the fundamental cosmological parameters which describe 
our observable universe (Pahcall et al. 1999| , [Eisenstein, Hu, fc Tegmark 1999| , P?urner 5: 



[lyson 1999| , |Wang, Spergel, fc Strauss 1999| ). The use of astrophysical standard candles 



provides a fundamental means of measuring the cosmological parameters Hq (current 
expansion rate of the universe), VLm (matter density fraction of the universe), and VL\ 
(vacuum energy density fraction of the universe). 

Type la supernovae (SNe la) are currently our best candidates for standard candles 
They can be calibrated to have small intrinsic dispersions ( [Phillips 1993| , [Riess, Press, fc 



[Kirshner 1995|) . Two independent teams, the High- 2; SN Search (Schmidt et al.) and the 
Supernova Cosmology Project (Perlmutter et al.), have observed about 100 SNe la. The 
data analysis results of both teams seem to suggest that our universe has a significant 
vacuum energy content ( [Garnavich et al. 1998| , [Riess et al. 1998| , Perlmutter et al. 1999|) . 



In this paper, we combine the data of the High- 2; SN Search team and the Supernova 
Cosmology Project, and bin the combined data by flux- averaging in redshift intervals. 
Previous work ( [Wang 2000| ) has shown that fiux-averaging is effective in reducing the 
scatter in SN la peak absolute luminosity due to weak gravitational lensing and intrinsic 
dispersions. Here, we study the effect of fiux-averaging on the estimation of cosmological 
parameters. In §2 we compare the data sets from the two teams. In §3 we describe how we 
combine the data from the two teams. In §4 we fiux-average the combined data and derive 
estimated cosmological parameters. In §5 we illustrate the effect of evolution of SN peak 
absolute luminosity on the estimation of cosmological parameters. §6 contains a summary. 



2. Comparison of data sets 



The published data of the High-z SN Search team ([Schmidt et al. 1998[ [Riess et"aL 



1998| ) consists of 50 SNe la. They give measured distance modulus for each SN la, /io, to 



be compared with the theoretical prediction 

where (^^(z) is the luminosity distance; it is related to the angular diameter distance dA{z) 
and comoving distance r{z): 

dLiz) = (1 + zfdAiz) = il + z) r{z), (2) 

assuming a completely smooth universe. 
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The comoving distance r{z) is given by ( [Weinberg 1972| ) 



riz) = 7—— -TTTr smn 



/ dz' 



-1/2 



(3) 



where "sinn" is defined as sinh if Qk > 0, and sin if Qk < 0. If Qk = 0, the sinn and QkS 
disappear from Eq.(0), leaving only the integral. 

The published data of the Supernova Cosmology Project ( Perlmutter et al. 1999| ) 
consists of 60 SNe la. They give the estimated effective B-band magnitude of each SN la, 
to be compared with 

= Mb + /ig, (4) 

where Mb is the peak absolute magnitude of a "standard" SN la in the B-band, and /Xg is 
given by Eq.(p. 

We find that the published data of the two teams have 18 SNe la in common. The two 
data sets for the same SNe la should differ by the constant offset Mb- We find 

M*-^^^^ = m^^^-/if^^^ = -19.33 ±0.25, (5) 
M]^i5 = m'^Jf - = -19.42 ± 0.27, (6) 

where /iQ^^*""^ and fi]^^^ for each SN la are estimated using the MLCS method and the 
template-fitting method (ml5) respectively by Riess et al. (1998). The difference in the 
data from the two teams has a dispersion of about 0.25 magnitudes. This is surprisingly 
large since 16 of these SNe la were drawn from the Hamuy et al. (1996) data. For these 16 
SNe la, we find m^^^ - /i*^^^^ = -19.30 ± 0.24, and m^^^ - /iJJ^^^ = -19.40 ± 0.26. For the 
rest of this section, we only consider these 16 SNe la. 

The simplest way to calibrate SNe la is to use the linear relation between maximum 
peak luminosity and decline time found by Phillips (1993). Hamuy et al. (1996) found 
M^„^. = a + b[Ami5{B) - 1.1], with a = -19.258(0.048) and b = 0.784(0.182) for 26 "low 
extinction" SNe la. Table 1 lists 16 SNe la from Hamuy et al. (1996) data that have been 
reanalyzed by both Riess et al. (1998) and Perlmutter et al. (1999). 

In Table 1, M^^^ is the corrected B band peak absolute magnitude given by Hamuy et 
al. (1996): 

ME'' = + 51og(iJo/65) - 0.784 {Arm, - l-l)- (7) 

We find that for the 16 SNe la listed, Mf = -19.253 ± 0.190. 
From Table 1, we find 

(m^/ - /if ^^^) - Mf 96 = -0.047 ±0.270 
(m^^^ - /iJT'^) - Mf 96 = -0.144 ±0.229. (8) 
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Table 1: 16 SNe la from Hamuy et al. (1996) 



SN la 


z 




— 


P f f 1 ^ 

rn^ — /ig 


A f r c n 

H'O ~ PO 


. . -1 pr Ti 

,,"il5 ,.P 

PO PO 


ini^ - fil 


1990af 


0.050 


-19.31 


-18.90 


-19.04 


-0.39 


-0.25 


-19.29 


1992P 


0.026 


-19.16 


-19.68 


-19.51 


0.31 


0.14 


-19.37 


1992ae 


0.075 


-19.21 


-19.37 


-19.34 


-0.05 


-0.08 


-19.42 


1992ag 


0.026 


-19.05 


-19.09 


-19.25 


-0.08 


0.08 


-19.17 


1992al 


0.014 


-19.48 


-19.45 


-19.66 


-0.16 


0.05 


-19.61 


iyy/aq 


U. iUi 


1 Q 1 7 


-iy .zo 


1 Q 1 7 

-iy. i ( 


n 1 /I 

-U. i4 




1 Q "iQ 

-iy.oy 


1992bc 


0.020 


-19.46 


-19.69 


-19.59 


0.00 


-0.10 


-19.69 


1992bg 


0.035 


-19.40 


-19.60 


-19.83 


0.14 


0.37 


-19.46 


1992bh 


0.045 


-18.85 


-19.30 


-19.26 


0.23 


0.19 


-19.07 


1992bl 


0.043 


-19.45 


-19.07 


-19.34 


-0.32 


-0.05 


-19.39 


1992bo 


0.018 


-19.22 


-19.11 


-19.27 


0.08 


0.24 


-19.03 


1992bp 


0.080 


-19.57 


-19.38 


-19.69 


-0.35 


-0.04 


-19.73 


1992br 


0.087 


-19.12 


-18.93 


-18.81 


0.01 


-0.11 


-18.92 


1992bs 


0.064 


-18.98 


-19.37 


-19.39 


0.12 


0.14 


-19.25 


19930 


0.052 


-19.32 


-19.49 


-19.77 


0.02 


0.30 


-19.47 


1993ag 


0.050 


-19.27 


-19.11 


-19.42 


-0.12 


0.19 


-19.23 



This indicates that the absolute magnitude of a "standard" SN la derived from comparing 
the Riess et al. (1998) and the Perlmutter et al. (1999) data sets differs appreciably from 
that derived by Hamuy et al. (1996) from the same 16 SNe la. 

To account for these differences, we examine the internal dispersion of each data set. 
Using fj,Q from Eq.(|l|) with z and Hq = 65kms~^ Mpc^^, we find 

/if^'^^-Zig = -0.043 ±0.195 
^l^^^-fil = 0.054 ±0.178 

rrf^S - = -19.34 ±0.22 (9) 

These should be compared with M^^^ = —19.253 ± 0.190 for the 16 SNe la calibrated 
using Eq.(^. Fig.l shows the internal dispersions in the calibrated SN la peak absolute 
magnitudes, as given by (a) Hq^^'^^ - f4, (b) fx]^^^ - f4, (c) m^J^ - /xg, and (d) Mf 
Clearly, the internal dispersion for each data set is about aint ~ 0.20 magnitude. This 
indicates that these 16 SNe la can be calibrated to be standard candles with a dispersion 
(intrinsic and observational) of about 0.2 magnitude. 

Next, we examine how the difference in analysis techniques introduces uncertainty. 
For the Riess et. al 1998 data analyzed using MLCS method and temp late- fitting method 
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respectively, we find 

= -0.097 ±0.163. (10) 
For the Perlmutter et al. (1999) data and the Hamuy et al. (1996) data, we find 

(m^^^-/ig)-Mf 9^ = -0.09 ±0.147. (11) 

It seems that the difference in analysis technique typically introduces a uncertainty of 
o'tech ~ 0.15 magnitude. 

Therefore, the difference between independently analyzed data sets could be as large 
as cr = af^^ ± (T^^ch ~ 0.25. The differences between the Riess et al. (1998) data and the 
Perlmutter et al. (1999) data stem mostly from the internal dispersion in the peak absolute 
magnitudes of SNe la estimated from each data set, with a substantial contribution from 
the difference in analysis techniques. 

Eq. (^) indicates that the SN la absolute magnitude estimated using rrf^^ — ^i^^^^^ is 
closer to that found by Hamuy et al. (1996). We will use the MLCS method data from 
Riess et al. (1998) for the rest of this paper. 



3. Combination of Data 

We combine the data from the two teams by adding 42 SNe la from the Supernova 
Cosmology Project to the High- 2; SN Search data (the MLCS method results), leaving out 
the 18 SNe la from the Supernova Cosmology Project which are already included in the 
High-z SN Search data. This yields a combined data set of 92 SNe la. In the combined 
data set, we convert rrf^^ of the the SNe la which have been taken from the Supernova 
Cosmology Project to /i^^^^^ by using Eq.(|) with Mb = -19.33. Mb = -19.33 is the 
mean difference in the data [roB^ — ft^^^'-^^) of the 18 SNe la published by both teams [see 
Eq.(^]. Fig. 2 shows the difference between m^-^ and fi^^^^^ for these 18 SNe la. 

Fig. 3 shows the magnitude-redshift plots of the combined data set of 92 SNe la. 
The solid points represent 50 SNe la from Schmidt et al. data. The circles represent 42 
additional SNe la from Perlmutter et al. The dashed line is the prediction of the open cold 
dark matter model (OCDM) with Qm = 0.2 and Q\ = 0. The dotted line is the prediction 
of the standard cold dark matter model (SCDM) with = 1 and = 0. The two thick 
solid lines are flat models with a cosmological constant (ACDM): {Qm, ^a) = (0.3, 0.7) 
and (fi„, ^a) = (0.2, 0.8). The thin solid line is the SCDM model (fi^, Ha) = (1, 0) with 
(1 ± z) dimming of SN la peak luminosity (linear evolution). Fig. 3(a) shows the distance 
modulus versus redshift. Fig. 3(b) shows the distance modulus relative to the prediction of 
the OCDM model. 
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3.1. Estimation of the cosmological parameters 

We find that the best-fit values of the cosmological parameters are sensitive to the 
estimated errors of the data points. Following Riess et al. (1998), we use 

r 1 2 

2/rr o ON /^0,i(^i 1^0, ^m, ^a) " /iO,i 
x\Ho,nm,^A) =Y.- 2 , 2 (12) 

where a^^ is the estimated measurement error of the distance modulus, and amz is the 
dispersion in the distance modulus due to the dispersion in galaxy redshift, cr^, due to 
peculiar velocities and uncertainty in the galaxy redshift (for the Perlmutter et al. data, 
the dispersion due to peculiar velocities is included in a eff, i.e., o"„„). Since 

B 

^-^=haoi^^j^- ^^^^ 

amz depends on Qa and Q^- We compute amz iteratively while estimating the cosmological 
parameters. 

For the Schmidt et al. data, we follow Riess et al. (1998) in adopting az = 200kms~^, 
and add 2500 kms~^ in quadrature to az for SNe la whose redshifts were determined from 
broad features in the SN spectrum. For the Perlmutter et al. data, we take a^^ = cr^eff, 

which already includes dispersion due to peculiar velocities of 300kms~^, and we use the 
redshift uncertainty az for each SN la given in their tables to compute amz (see Eq.([T^)). 

Table 2 lists the derived best-fit cosmological parameters, with 1-a error bars, for the 
High-z SN Search data (Schmidt et al.), the Supernova Cosmology Project data (Perlmutter 
et al.), and the combined data. Note that h is the dimensionless Hubble constant, defined 
by Hq = 100 /ikms~^Mpc~^. xl P^r degree of freedom. 

Fig.4 shows the 68.3% and 95.4% confidence contours in the — Qm plane. The 
dotted lines represent the Schmidt et al. data (50 SNe la), the dashed lines represent the 
Perlmutter et al. data (60 SNe la), and the solid lines represent the combined data (92 SNe 
la). Note that we have allowed Qm to be a free parameter in estimating its value. Since 
distances are not directly measurable, we can treat them as theoretical intermediaries in 
the data analysis, therefore Qm must be allowed to have negative as well as positive and 
zero values, for statistical robustness. 

We have combined data from two independent teams analyzed using different statistical 
methods. It is important to analyze all the SN data with the same techniques for statistical 
consistency; this is not possible at present, because not all the reduced SN la data have 
become public (Perlmutter et al. have not yet published the light curves of their 42 SNe la). 
However, the published data from the two teams are similar enough (with a constant offset) 
to allow a meaningful study of the effect of flux- averaging on the estimation of cosmological 
parameters from SN la data (see §4 and §5). 
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Table 2: Estimated cosmological parameters with SNe la 





50 SNe la (Schmidt et al.) 


60 SNe la (Perlmutter et al.) 


92 SNe la (combined data) 


h 


0.65 ± o.or 


0.66 ± 0.02* 


0.65 ± 0.01* 




0.2 ± 0.6 


1.0 ± 0.4 


0.7 ± 0.4 


^2 


0.7 ± 0.8 


1.7 ± 0.6 


1.2 ± 0.5 


"X.V 


1.16 


1.63 


1.48 



Note. — * Statistical error only, not including the contribution from the much larger SN 
la absolute magnitude error. 

3.2. Hq as a systematic indicator of data sets 

We have chosen to estimate Hq simultaneously with VLm and mainly for three 
reasons. 

First, Hq is a useful indicator of the systematic "zero point" of a given SN la data 
set. Let us consider the 16 SNe la from Hamuy et al. (1996) that both teams (Riess et 
al. 1998, Perlmutter et al. 1999) have reanalyzed (see §2). Since each of these 16 SNe la 
can only have one true peak absolute mag nitude, we take M^^^^ = = Mf where 

M^LCS^ j^ejf ^ j^H96 ^^le peak absolute magnitudes derived from the Riess et al. 
1998 (MLCS), Perlmutter et al. 1999, and Hamuy et al. 1996 data respectively. From 
Eq.(il), /i^^^^^ - /io = -0.043 and m^'J^ - f-i^ = -19.34. Using m'J^ = M^J^ + fx'J^ , and 
setting M^J^ = Mf = -19.253, we find fi^^^ - /io = -0.087. Hence 

/io*^"^^-/^o^^ = -1^^ = 0.044, (14) 

i.e., the difference between the estimated values for the Hubble constant from Riess et al. 
1998 {H^^^^^) and Perlmutter et al. 1999 {Hq^^) data is 

^MLC5_^e//^ _^_32^ (15) 

for Hq^^^^ = 65. This is in agreement with the estimated values of Hq listed in Table 2 for 
50 SNe la from the Riess et al. (1998) data and for 60 SNe la from the Perlmutter et al. 
(1999) data. The SNe la from Perlmutter et al. (1999) are systematically brighter than the 
Riess et al. (1998) SNe la by about 0.04 magnitude. 

Second, integrating out the Hq dependence in the probability distribution functions of 
the cosmological parameters have very little effect on the estimated values of Qm and f^A, 
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because the dependence is uncorrelated with the and f^A dependence. Therefore, 
we have httle to gain in the accuracy of estimated Qm and by excluding Hq from the 
parameter estimation. 

Thrid, including Ho in the parameter estimation conforms with the common practice 
of estimating all basic cosmological parameters simultaneously when analyzing other 
cosmological data sets, for example, the Cosmic Microwave Background Anisotropy data 
( |Wang, Spergel, fc Strauss 1999| , pHisenstein, Hu, fc Tegmark 1999| ). The distance modulus 



depends on Hq, as well as f2m and f^A. Therefore, Hq should be estimated from the 
SNe la data simultaneously with and f^A, although the Hq dependence of the data is 
independent of the dependence on Qm and f^A- 

To summarize, including Hq in the cosmological parameter estimation from SN la data 
provides a useful systematic indicator of the data set used, while having little effect on 
the accuracy of the other cosmological parameters [Qm and ^a) being estimated, and it 
conforms with the convention used in parameter estimation from other cosmological data. 
Note that the error of the Hq estimated here only reflects statistical error, it does not 
include the contribution from the much larger SN la peak absolute magnitude error. For a 
realistic assessment of the actual value and errors associated with Hq estimated from SNe 
la, the reader is referred to Saha et al. (1999), Jha et al. (1999), and Gibson et al. (1999). 



4. Flux-averaging of data 



An important reason to consider flux- averaging is the weak gravitational lensing of 
SNe la. Because our universe is inhomogeneous in matter distribution, weak gravitational 
lensing leads to a non-Gaussian distribution in the magnification of standard candles. 

For a given redshift z, if a mass-fraction a of the matter in the universe is smoothly 
distributed, the largest possible distance for hght bundles which have not passed through a 
caustic is given by the solution to the Dyer-Roeder equation ( Pyer fc Roeder 1973| , Schneider 
et al. 1992|, ^antowski 19980: 



d 

9[z) — 
dz 



dD, 



dz 



Da{z = 0) = 0, 



+ -aQU'^ + zfDA = 0, 



dD, 



dz 



c 

H'q 



(16) 



where g{z) = (1 + z)'^y 1 + VlmZ + Vt\[{l + z)^'^ — 1]. The smoothness parameter a 
essentially represents the amount of matter that can cause the magnification of a given 
source. If we define a direction-dependent smoothness parameter a via the Dyer-Roeder 
equation, there is a unique mapping between a and the magnification of a source. We 
can think of our universe as a mosaic of cones centered on the observer, each with a 
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different value of a. Wang (1999) has derived empirical fitting formulae for the probability 
distribution of a from the numerical simulation results of Wambsganss et al. (1997). 

Wang (2000) showed that flux- averaging of simulated data (with noise due to 
both weak lensing and intrinsic dispersions) leads to SN peak luminosities which well 
approximate the true luminosities with a = \ (which represents a completely smooth 
universe). The angular-diameter distance defined in Eq.(|^) using the comoving distance 
r{z), dA^z) = r{z)/{l + z), satisfies Eq.([l6|) with (5 = 1, i.e., the angular-diameter distance 
defined in Eq.(0) is the Dyer-Roeder distance with a = 1. Therefore, flux- averaging justifies 
the use of Eq. (0) in the analysis of SN la data. 

Before fiux-averaging, we convert the distance modulus of SNe la into "fluxes", 
f{zi) = l0-^^oizi)/2.5 _ then obtain "absolute luminosities", {C{z,i)}, by removing the 
redshift dependence of the "fluxes", i.e.. 



where is the luminosity distance, and {HQ,flm,^A) are the best-fit cosmological 
parameters derived from the unbinned data set {f{zi)}. We then flux-average over 
the "absolute luminosities" {Ci} in each redshift bin. The set of best-fit cosmological 
parameters derived from the binned data is applied to the unbinned data {f{zi)} to obtain 
a new set of "absolute luminosities" {Ci}, which is then flux-averaged in each redshift bin, 
and the new binned data is used to derive a new set of best-fit cosmological parameters. 
This procedure is repeated until convergence is achieved. The l-cr error on each binned data 
point is taken to be the root mean square of the 1-a errors on the unbinned data points in 
the given redshift bin, {/j} {i = 1,2, A^), multiplied by (see [Wang 20^0]) . 



Fig.5 shows magnitude-redshift plots of the binned data for the total of 92 SNe la, 
with redshift bin Az = 0.05. The lines are the same as in Fig.l. Fig. 5(a) shows the distance 
modulus versus redshift. Fig. 5(b) shows the distance modulus relative to the prediction of 
the open cold dark matter model (OCDM) with = 0.2 and Q/^ = 0. Fig. 6 is the same as 
Fig.5, but with redshift bin Az = 0.1. 

Table 3 lists the estimated cosmological parameters, with 1-a error bars. 

Fig. 7 shows the 68.3% and 95.4% confidence contours in the Q\ — Qm plane. The solid 
lines represent the unbinned data, the dotted lines represent the binned data with redshift 
bin Az = 0.05; the dashed lines represent the binned data with redshift bin Az = 0.1. 

Fig. 8 shows that the estimated parameters from flux- averaged data are not sensitive to 
the size of the redshift bin. The thick solid lines are the estimated flm (a) and (b) as 
functions of the redshift bin size dz, the dotted lines mark the la errors on (a) and Q\ 
(b). The thin solid line is xl + 2. The wiggles in the lines are due to the small number of 
SNe la in each redshift bin for small bin size, and due to the small number of binned data 
points for large bin size. Clearly, the optimal range for the bin size is 0.025 < dz < 0.9. 



C{Zi) = 4:7ldl{Zi\Ho,nm,^A) f{Zi) 



(17) 
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Table 3: Estimated cosmological parameters with 92 SNe la 





unbinned data 


binned with Az — 0.05 


binned with Az — 0.1 


N, . 


92 


17 


10 


h 

n 


0.65 ± 0.01* 


0.65 ± 0.01* 


0.65 ± 0.01* 


r) 


0.7 ± 0.4 


0.3 ± 0.6 


0.3 ± 0.6 




1.2 ± 0.5 


0.7 ± 0.7 


0.7 ± 0.7 


xl 


1.48 


0.73 


0.78 




fixing Q\ = 1 — 


I. 
tl 


0.65 ± o.or 


0.65 ± 0.01* 


0.65 ± 0.01* 




0.3 ± 0.1 


0.3 ± 0.1 


0.3 ± 0.1 




1.47 


0.68 


0.68 




fixing Q\ = 


h 


0.64 ± 0.01* 


0.65 ± 0.01* 


0.65 ± 0.01* 




-0.2 ± 0.1 


-0.2 ± 0.1 


-0.2 ± 0.1 


x' 


1.49 


0.72 


0.76 




fixing flm = 0.3 


h 


0.65 ± 0.01* 


0.65 ± 0.01* 


0.65 ± 0.01* 




0.7 ± 0.1 


0.7 ± 0.2 


0.7 ± 0.2 


x' 


1.47 


0.68 


0.68 



Note. — * Statistical error only, not including the contribution from the much larger SN 
la absolute magnitude error. 

Because fiux-avcraging reduces the bias due to weak lensing, the flux-averaged data 
yield less biased estimates of the cosmological parameters. 



5. Effect of evolution 

Drell, Lorcdo, & Wasserman (1999) have studied the effect of SN la luminosity 
evolution on the estimation of cosmological parameters. Recently, Riess et al. (1999) 
suggested that there is an indication of evolution of SNe la from their risetimes, and 
that if this observed evolution affects the peak luminosities of SNe la, it can account for 
the observed faintness of high- 2; SNe la without invoking a cosmological constant. To 
illustrate this interesting possibility, we introduce a power-law evolution in the SNe la peak 
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luminosity, C{z) = Ati d\{z) f{z), i.e., 

C{z) = {l + z)'^C{z = 0). (18) 

Here, /3 > represents brightening of SN la with redshift, while f3 < represents dimming 
with redshift. The thin solid line in Fig. 3, Fig.5, and Fig.6 represent the SCDM model 
{Qrn = 1; = 0) wlth linear dimming of the SN la peak luminosity with f] = —1. 

Table 4 lists the estimated cosmological parameters and SN la evolution power-law 
index f3, with 1-a error bars, for the combined data set of 92 SNe la. 

Table 4: Estimated cosmological parameters with SN evolution 





unbinned data 


binned with Az = 0.05 


binned with Az = 0.1 


^data 


92 


17 


10 


h 


0.66 ± 0.02* 


0.65 ± 0.02* 


0.65 ± 0.03* 




1.5 ± 5.8 


0.8 ± 22.5 


0.9 ± 31.3 




1.2 ± 0.8 


0.4 ± 4.5 


0.4 ± 3.3 


/? 


-0.7 ± 3.7 


-0.6 ± 17.1 


-0.6 ± 21.1 




1.49 


0.78 


0.91 




fixing Qm = 0.3 


h 


0.65 ± 0.01* 


0.65 ± 0.01* 


0.65 ± 0.01* 




1.4 ± 1.1 


0.7 ± 1.3 


0.6 ± 1.4 




0.5 ± 1.6 


0.0 ± 1.0 


0.0 ± 1.0 




1.59 


0.73 


0.78 



Note. — *Statistical error only, not including the contribution from the much larger SN 
la absolute magnitude error. 

There is strong degeneracy, as expected, between the evolution power-law index (3 
and the matter density fraction Qm- If the SNe la peak absolute luminosity changes with 
redshift due to evolution, our ability to measure the cosmological parameters from SN la 
data will be significantly diminished, unless we can correct for the evolution. The issue of 
SN evolution can only be settled if a large number of SNe la are observed a.t z > 1; this can 
be accomplished via a supernova pencil beam survey using a dedicated 4 meter telescope 
( IWang 200q ). 
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6. Conclusions 

We have combined the data from the two independent SN la groups, Schmidt et al. 
and Perlmutter et al., and analyzed the combined data using flux-averaging in redshift bins 
of 0.05 and 0.1. We flnd that the estimation of the cosmological parameters is not sensitive 
to the size of the redshift bin. While the combined data without flux-averaging are best 
flt by a closed universe with high matter content and larger than critical density vacuum 
energy, the flux-averaged data are best flt by a nearly flat universe with a low matter 
content. This is consistent with the strong observational evidence that we live in a low 
matter density universe ([Bahcall, Lubin, fc Dorman 1995| , [Carlberg et al. 1996| , Pahcall ^ 



Fan 1998] , [Krauss 1998| , [Bahcall et al. 19^9] ). Flux-averaging leads to less biased estimates 



of the cosmological parameters by reducing the bias due to systematic effects such as weak 
lensing. 

The distance- redshift relation, Eq.(3), is generally used in making theoretical 
predictions to compare with the brightnesses of standard candles. However, it is only 
valid in a smooth universe. We live in a clumpy universe. Because of flux conservation, 
flux-averaging justifles the use of Eq.(3) in the analysis of SN la data (see §4). 

In the estimation of cosmological parameters, we have allowed the matter density 
fraction and the vacuum energy density fraction Q\ to be unconstrained free parameters. 
Unlike fluxes, distances are not directly measurable, therefore they are best treated as 
theoretical intermediaries. Since the cosmological parameters Qm and Qa enter through 
distances [see Eqs.(2) & (3)], they should be allowed to have negative as well as positive 
and zero values for statistical robustness. By not applying priors on and Q\ in our 
parameter estimation, we can check the validity of Eqs.(2) & (3). For example, if the most 
probable value of Qm derived from the data is negative {Q^ < 0), it could be an indication 
that Eqs.(2) and (3) do not apply, since we know that Q^n > 0. 

We have also investigated the effect of possible evolution of SNe la on the estimation 
of cosmological parameters. Assuming power-law evolution in the SN la peak absolute 
luminosity, we flnd that there is strong degeneracy between the evolution power-law index 
and the matter density in the universe. The evolution of SNe la must be resolved before 
they can be regarded as reliable standard candles. Supernova pencil beam surveys which 
could yield hundreds of SNe la at z > 1 over a few years will be critical in constraining the 
evolution of SNe la ( |Wang 2000| ). 
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Fig. 1. — The internal dispersions in the cahbrated SN la peak absolute magnitudes from 
four data sets, (a) //^^^^ - //g, (b) //^^^ _ (^-) ^ejf _ ^Hm_ 

Fig. 2. — The difference between m^J^ (Perlmutter et al 1999) and /i^^'^^ (Riess et al 
1998) for the same 18 SNe la. The error bars are the combined errors in mfj-^ and fj,^^^^. 

Fig. 3. — The magnitude-redshift plots of the combined data set of 92 SNe la. The solid 
points represent 50 SNe la from Schmidt et al. data. The empty points represent 42 
additional SNe la from Perlmutter et al. data. The dashed line is the prediction of OCDM 
with Q„i = 0.2 and Q\ = 0. The dotted line is the prediction of SCDM with Qm = 1 and 
Qa = 0. The two thick solid lines are predictions of ACDM, with {Qm, ^a) = (0.3, 0.7), 
{flm, ^a) = (0.2, 0.8) respectively. The thin solid line is SCDM with (1 + z) dimming of 
SN la peak luminosity, (a) The distance modulus versus redshift. (b) The distance modulus 
relative to the prediction of the OCDM model. 

Fig. 4. — The 68.3% and 95.4% confidence contours in the Q\ — Qrn plane. The dotted lines 
represent the Schmidt et al. data (50 SNe la), the dashed fines represent the Perlmutter et 
al. data (60 SNe la), and the sofid lines represent the combined data (92 SNe la). 

Fig. 5. — The magnitude-redshift plots of the binned data for 92 SNe la, with redshift bin 
Az — 0.05. The fine types are the same as in Fig.3. (a) The distance modulus versus 
redshift. (b) The distance modulus relative to the prediction of the open cold dark matter 
model (OCDM) with = 0.2 and Qa = 0. 

Fig. 6. — The same as Fig.5, but with redshift bin Az — 0.1. 

Fig. 7. — The 68.3% and 95.4% confidence contours in the Qa — plane. The solid lines 
represent the unbinned data, the dotted lines represent the binned data with redshift bin 
Az — 0.05; the dashed lines represent the binned data with redshift bin Az = 0.1. 

Fig. 8. — The estimated parameters from fiux-averaged data as functions of the size of the 
redshift bin. The thick sofid fine is the estimated parameter, with la errors marked by the 
dotted lines. The thin solid fine is + 2. (a) Q^, (b) 0,a- 
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Fig. 1. — The internal dispersions in the cahbrated SN la peak absolute magnitudes from 
four data sets, (a) n^^'^^ - n^. 
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Fig. 1.- (b) _ 
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Fig. 1.— (c) m'^J^ - /xg. 
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Fig. 1.— (d) Mf 
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Fig. 2. — The difference between rrf^^ (Perlmutter et al. 1999) and Hq^^^'^ (Riess et al. 
1998) for the same 18 SNe la. The error bars are the combined errors in m'^J^ and n^^^^ . 
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Fig. 3. — The magnitude-redshift plots of the combined data set of 92 SNe la. The sohd 
points represent 50 SNe la from Schmidt et aL data. The empty points represent 42 
additional SNe la from Perlmutter et al. data. The dashed line is the prediction of OCDM 
with Qrn = 0.2 and fl\ = 0. The dotted line is the prediction of SCDM with = 1 and 
Q\ = 0. The two thick solid lines are predictions of ACDM, with {Qm, ^a) = (0.3, 0.7), 
{flm, ^a) = (0.2, 0.8) respectively. The thin solid line is SCDM with (1 + z) dimming of SN 
la peak luminosity, (a) The distance modulus versus redshift. 
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Fig. 3. — (b) The distance modulus relative to the prediction of the OCDM model. 
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Fig. 4. — The 68.3% and 95.4% confidence contours in tlie fl\ — flm plane. The dotted lines 
represent the Schmidt et al. data (50 SNe la), the dashed hnes represent the Perlmutter et 
al. data (60 SNe la), and the sohd lines represent the combined data (92 SNe la). 
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Fig. 5. — The magnitude-redshift plots of the binned data for 92 SNe la, with redshift bin 
= 0.05. The hne types are the same as in Fig.3. (a) The distance modulus versus 
redshift. 
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Fig. 5. — (b) The distance modulus relative to the prediction of the open cold dark matter 
model (OCDM) with = 0.2 and Qa = 0. 
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Fig. 6. — The same as Fig. 5, but with redshift bin iS.z = 0.1. (a) The distance modulus 
versus redshift. 
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Fig. 6. — (b) The distance modulus relative to the prediction of the open cold dark matter 
model (OCDM) with = 0.2 and Qa = 0. 
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Fig. 7. — The 68.3% and 95.4% confidence contours in tlie fl\ — flm plane. Tlie solid lines 
represent the unbinned data, the dotted lines represent the binned data with redshift bin 
Az — 0.05; the dashed lines represent the binned data with redshift bin Az = 0.1. 
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Fig. 8. — The estimated parameters from flux-averaged data as functions of the size of the 
redshift bin. The thick sohd hne is the estimated parameter, with la errors marked by the 
dotted hnes. The thin sohd hne is + 2. (a) Q^. 
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Fig. 8.- (b) Qa. 



